Eigenvalues and Eigenvectors of the Staggered Dirac Operator at Finite Temperature 
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We examine the eigenvalues and eigenvectors of the staggered Dirac operator on thermal ensembles 
created in QCD with two flavours of staggered quarks. We see that across the phase transition a 
gap opens in the spectrum. For finite volume lattices in the low-temperature phase the eigenvectors 
are extended, but generic field configurations in the high temperature phase give rise to localized 
eigenstates. We examine measures of the stability of such localization and find that at finite volumes 
localization occurs through Mott's mechanism of the formation of mobility edges. However, the band 
gap between the localized and extended states seem to scale to zero in the limit of large volume. 
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I. INTRODUCTION 



Any fermionic operator can be written in the spectral form 



j§: o = j20ux\X){»\, (i) 

where | A) is an eigenvector of the Dirac operator with eigenvalue A, evaluated separately on each configuration. Typical 
operators of interest contain quark loops with various insertions, i. e., O = Tr {A\{D + m)~ 1 A2(D + m) -1 ■ ■ ■ A n {D + 
to) -1 ). As a result, 

n _ ST (Ai|^i|A2)(A 2 |A 2 |A 3 )---(A n |A 2 |A 1 ) 

o 

where we use the symbol D to refer to the massless Dirac operator. If the Ai commute with D, then the matrix 
elements in the numerator are diagonal, and all questions about the operator reduce to the simultaneous eigenvalues 
of the Ai and the Dirac operator. This happens, for example, in the chiral sector of the theory, where one deals 
• i-H , with questions about n-point functions of pions. Since 75-D75 = D\ most questions about the chiral sector can be 
answered if the eigenvalues are known. As a result, the thrust of many previous studies of QCD to date has been on 
the spectrum of eigenvalues, particularly on comparisons with random matrix theory (RMT) This focus is due to 
the fact that RMT is known to be equivalent to chiral perturbation theory in some limits 

However, at finite temperature, especially above T c , chiral perturbation theory is not the appropriate long-distance 
effective theory. Furthermore, there are interesting questions at many different length scales and one may need to 
build different effective theories to answer these questions. Several questions involve fermionic loops with insertions 
of operators which do not commute with D. An example is the vector susceptibility, 

v _ v K A ^|A 2 )I 2 f - 

Xv 2^ (m + AiKm + Aa)' W 

which includes quark number susceptibilities. Deeper understanding of such quantities need the study of the eigen- 
vectors 0. 
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II. EIGENVALUES 



We analyzed configurations generated in the study of QCD with two flavours of dynamical staggered quarks at a 
lattice spacing a = 1/4T |. The scale fixing yielded T c /m p = 0.186 ± 0.006. As T varied between 0.75T C and 2T C , 
the renormalized quark mass was kept constant. The physical box size, L = N s a where N s is the box size in units of 
the lattice spacing. The aspect ratio was varied in the range 2 < LT < 6. 

We investigated the eigenvalues, A, and eigenvectors, |A), of the massless Dirac operator, D, in typical thermal 
ensembles picked from these simulations. We used five configurations separated by two autocorrelation times at all 
temperatures and volumes except at 1.05T C where we verified the results using twenty configurations. Eigenvalues 
and eigenvectors were computed with the ARPACK subroutines [5]. For convergence, the tolerance is chosen so that 

|r| 2 < e, where r = (D — (4) 

where A is an eigenvalue and tp\ is the corresponding eigenvector. We report results with e — 2 x 10~ 13 . 
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FIG. 1: The lowest Dirac eigenvalues as a function of T for different spatial lattices. In all cases the lattice spacing a = 1/4T. 
Also shown in the lowest Matsubara frequency expected at this cutoff, i. e., the expectation for free fermions. 

The lowest staggered Dirac eigenvalue, Ao, evolves as shown in Figure [T] There is a clear crossover from low to 
high temperature behaviour evidenced by an increase in the lowest eigenvalue by three order of magnitude in the 
neighbourhood of T c . This becomes sharper in the neighbourhood of T c with increasing spatial size of the lattice. 

It is also worth noting that at 2T C there seems to be little remaining volume dependence. It is interesting that this 
large volume behaviour sets in at a minimum L given by LXq w 0.8. Consistent with this observation, at 1.5T C the 
lattice sizes which satisfy this condition also give results which are almost volume independent. 

In Figure [T] we further show the lowest eigenvalue of the free Dirac operator on a lattice of the same coarseness. 
In the limit of zero lattice spacing this would correspond to the Matsubara frequency, = irT. The full QCD 
configurations can be seen to lie very far from the free field theory (ideal gas) limit at all temperatures up to 2T C . 
This is consistent with another observation at the same lattice spacing that the pseudoscalar screening correlator 
constructed with staggered Dirac quarks yields a screening mass far below that expected from the free theory @. It 
would be interesting to see whether this correlation, related to the chiral behaviour, changes in the same way due to 
various improvements in the gauge and fermion actions and in the continuum limit. 

The Banks-Casher formula Q relates the density of Dirac eigenvalues at zero with the chiral condensate. In order 
to utilize this formula we expand the cumulative distribution of the eigenvalues in the form 

I(x) = / d\p(X) = a n x n , from which p(x) = na n x n ~ l . (5) 

^° n>l n>l 
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TABLE I: Fits of the cumulative density for (3 = 5.26, i. e., T/T c = 0.90 ± 0.01. 

where the 75-Hermiticity of the staggered Dirac operator as well as its anti-Hermitean nature have been used. Together 
they imply that the eigenvalues are paired and imaginary, The integration above is over the positive A values. 

Note that the reflection symmetry in A permits the existence of an a-i term (even n, in general) only if the p(A) is 
non-analytic at the origin. 

The cumulative distribution was constructed numerically and fitted to the form in eq. (JSJ) for configurations above 
and below T c . Indicative results are shown in Tables III and ITU The tables are arranged in increasing order of N, the 
number of eigenvalues included, and n of eq.©. Note that below T c one gets a good determination of a\ for cutoffs 
of the order of 0.1 or so. In fact, the values of ai do not depend on n or N. With increasing lattice size, L, the 
estimate of a\ increases marginally. On the larger lattices 02 is compatible with zero, indicating that the distribution 
is analytic. 

Above T c , 0,2 is clearly non-zero for all cut-offs while a± drops with increasing L. This behaviour, shown in Table 
Inl implies that a non-analyticity develops in the spectral density. This non-analyticity is due to the formation of a 
gap — the spectral density is exactly zero upto the gap, and then becomes non-zero. Another way to test this would 
be to introduce a gap explicitly in the eq. ([5]), 

p(x) = ^na n (x - xoY 1 ^ 1 . (6) 

n>l 

Indeed, when one does that, a non-zero value of the gap, xq, is observed for those temperatures where ai is non-zero 
by the other method. 

III. EIGENVECTORS AND MEASURES OF LOCALIZATION 

The eigenvectors of the Dirac equation, ip, are often investigated through the localized moments 

^(A) = ^"- 1 ^b 7 (r;A)|", where p 7 (r; A) = (A| 7 |A), and P* = 1, (7) 

r 

7 is a matrix in Dirac space [13], the inner product in the definition of p 7 (r; A) involves a sum over spin-flavour and 
colour indices, the explicit sum is over all V lattice sites r, and the normalization of the eigenvectors, involves the 
density where the Dirac matrix is identity. For staggered quarks one has the identity P* = P^ 5 . The second moment, 
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TABLE II: Fits of the cumulative density for /3 = 5.30, i. e., T/T c = 1.05 ± 0.01. 



n = 2 is called the inverse participation ratio (IPR). The moments have the interesting property that for constant 
ip = 1/VV, one finds P^ = 1, whereas for the localized tp(r) = S r ^ rg , one has P* — V n_1 . 

Histograms of IPR against A are shown in Figure [51 There is a very clear difference between the IPR observed 
below and above T c . Below T c the IPRs are close to unity, without any clear dependence on the eigenvalues. In 
contrast, the situation is dramatically different above T c ; several eigenvectors have very large values of IPR. There is 
correlation between the eigenvalue and IPR, with larger eigenvalues coming with substantially smaller IPR. 

Below T c there is little sign of volume dependence of the IPR, consistent with the small values seen there. Above 
T c the smaller IPR values seen for large A are also volume independent. However, as shown in Figure^ larger values 
of IPR are volume dependent. A test of scaling shows that the lattice size dependence is consistent with a power 
behaviour, L a , with 2.5 < a < 3.5. Again, this is not unexpected, since IPR is constructed to be proportional to the 
volume for localized eigenvectors. 

In [3j the transition from volume dependent to independent values of IPR is used to locate the "mobility edge" . By 
this identification one would have a mobility edge at A ~ 1.25T C for a temperature of 2T C . However, the notion of a 
mobility edge contains more physics and we shall examine it more critically in a later section. 

The eigenvalues and eigenvectors of the Dirac operator are clearly dependent on the gauge field backgrounds. 
However, thermodynamic quantities constructed from these have fluctuations which decrease rapidly with increasing 
lattice size. The IPR is not such a variable: its fluctuations are comparable to the average, as can be seen in Figure [3] 
The ratio of the variance and mean of P£, as a function of A at 2T C , is of order unity [ll|. The localization properties 
of Dirac eigenfunctions can therefore serve to classify the ensemble of gauge configurations which give important 
contributions to the thermal path integral. This is an obvious statement for overlap quarks, where localized chiral 
eigenvectors of the overlap Dirac operator are closely connected to localized gauge field configurations which are taken 
to be the lattice analogue of instantons. It is interesting that localization using staggered quarks, where the connection 
to topology is obscure, can also be used as a tool for analysis of gauge configurations. 

The notion of localization has been closely examined in |8| . Since pi (r) is non-negative and normalized to unity 
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FIG. 2: The IPR, P 2 \ above and below T c as a function of the staggered Dirac eigenvalue A. 



one can construct a measure of localization in the following way. Take a value p / and find the fraction of the lattice 
sites, f(j>f), containing values pi(r) > Pf. Clearly f(pf) lies between and 1, and is a decreasing function of pf. The 
integral of Pi(r) over these sites, C(pf), lies between and 1, and is another decreasing function of pf. Eliminating 
Pf between these two, one obtains Horvath's localization function f(C). Clearly /(C = 0) = 0, f(C = 1) = 1 and the 
function is non-decreasing. 

If Px(f) is highly peaked, then C(pf) increases rapidly as pf decreases, whereas /(p/) increases slowly. As a result, 
/(C) is small over most of the range of C as increases very rapidly to unity near the end of the range. If, on the 
other hand, pi(r) is fairly uniform, then both C(pf) and f(pf) increase fairly abruptly over a small range of pf. The 
function /(C) then increases very rapidly towards unity at small C. In Figure [4] we show the behaviour of several 
models ofpi(r) — 
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FIG. 3: The relative fluctuations in the IPR, i. e., the ratio of the variance and the mean, at 2T C as a function of the staggered 
Dirac eigenvalue, A. 
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FIG. 4: Examples of the localization function, /(C), for different models of localization, compared with two examples from 
QCD. 



1. Some periodic functions, cos 2 (fc • r), normalized to unity on two lattices; these have P% = 1-5- 

2. Random uncorrelated function values on sites, drawn from the uniform distribution in [0, 1], normalized to unity; 
these have P£ = 1.66. 

3. Two Dirac eigenvectors obtained from the same gauge configuration at 2T C , one with P^ < 2 and the other with 
Pi > 100. 

The localization function f(C) clearly contains more information than the single number P£, i. e., the IPR. However, 
the IPR is statistically compatible with statements obtained from the more detailed measurement of f(C). We 
demonstrate this by the following correspondence. Choose any arbitrary value, C*, the function value /(C*) is strongly 
correlated with the IPR, as we show in Figure [5] For a wide range of C» we find /(C*) oc I/-P2 ■ 

We give an example of a question which can be easily answered through the use of the localization function. If the 
eigenvector is localized, then how does it fall off away from the peak? Exponential fall, ip(r) ~ exp(— aR), where R 
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FIG. 5: The strong correlation between IPR, P£, and the value of the localization function for two values of C shows that the 
latter contains all the information available in the former. 
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FIG. 6: Scaling of the localization function in the vicinity of / = C = 1 shows that staggered Dirac eigenvectors with IPR 
larger than 40 fall exponentially far from the peak, whereas those with IPR less than 2 have drastically different behaviour. 



is the distance from the peak, would imply C ~ 1 — g{R) exp(— R 2 ) and / ~ R d , in d dimensions. Thus, exponential 
falloff of a localized eigenvector would give rise to the relation — log(l — C) oc %/7i f° r both C and / close to unity. 
In Figure El we show that this is true of staggered Dirac eigenvectors with > 40 but those with P 2 X < 2 have 
completely different behaviour. 

A model for eigenvectors with small IPR is that of a function pi(r) with random uncorrelated values. We call this 
the site percolation model for the following reason. As we trace out the level curves of this function by choosing pf, 
we pick sites independently with a probability given exactly by /. Each site belongs to an unique cluster, defined as 
the collection of all neighbouring sites on the lattice which are picked [l2j ■ When / is small, we find small localized 
clusters, but above some critical value, we have percolating clusters. Each realization of the random function is a 
realization of the percolation problem for all possible probabilities. 

As we fill a larger and larger fraction of the lattice, the number of clusters, N c , grows until the percolation threshold 
is reached, after which the number of clusters begins to decrease. The clusters are ramified, and, near the critical 
percolation probability, have a fractal dimension related to the critical indices of the percolation problem. Above the 
critical porbability, the clusters have canonical dimension, as a result of which the holes are filled in rapidly, and N c 
decreases. 

In Figure [7] we compare the average cluster size as a function of / for the site percolation problem and those 
eigenvectors in QCD at 2T C which have IPR greater than 2. The fact that QCD has more clusters at larger / than 
site percolation implies that the percolating cluster constructed from P^ have larger holes inside them where isolated 
clusters can exist. 
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FIG. 7: The number of clusters, N c , normalized by the lattice volume for uncorrelated site percolation and for those eigenvectors 
of the staggered Dirac operator at 2T C which have IPR smaller than 2. 



IV. STABILITY OF LOCALIZATION 




FIG. 8: The stability of localization of staggered Dirac eigenvectors with 2 flavours of dynamical quarks in QCD at 2T C . The 
first panel shows that stability decreases with increasing spatial size. The second shows that the data supports scaling as 1 /L 2 . 

One of the paradigms in the analysis of Dirac eigenvectors is that of Mott localization and the existence of a 
mobility edge. In a metallic crystal with random impurities, localization of electron wavefunctions can be observed. 
Mott argued that if there exist a localized and an extended state arbitrarily close in energy, then they will mix under 
any small perturbation of the Hamiltonian (induced, for example, by the movement of one of the impurities) hence 
destroying localization. He argued that, as a result, localization is robust only when localized and extended states 
are separated in energy. It is well-known that this argument could fail if the extended states have support in regions 
with holes, since the lack of overlap can then be arranged in space rather than in energy. 

The mobility edges are the band edges of localized states. On a finite lattice where the eigenvalue spectrum is 
discrete, the identification of a mobility edge is not straightforward. As a result, it is hard to test Mott's picture of 
localization directly. It is interesting to build another measure of stability. We do this next. 

Assume that the Dirac operator is perturbed by a change in the gauge fields, D(U + 8U) = D(U) + 8D. Then, first 
order perturbation theory tells us that the change in an eigenvector is 

8\\)=Y,Cxv\\'), where C AA , = A > . (8) 
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FIG. 9: Correlation between the measures of stability of the localization of staggered Dirac eigenvectors, with 2 flavours of 
dynamical quarks in QCD at 2T C . Note the common scaling of the stability measure 5a with the eigenvalue difference Q\. 



Under a random change of the gauge fields, the phase information in the matrix element above is randomized. Hence, 
for a study of average properties of the perturbations over ensembles of random changes of gauge field, it would suffice 
to study C8U instead, where 



C 



|A- A*l 



(9) 



This matrix can be extracted purely from the knowledge of the eigenvalues and eigenvectors of the staggered Dirac 
operator. Note also that the mixing involves both a spatial part, which is the numerator, and a part in energy, which 
is the denominator. A small mixing can be a result of either. 

A perturbing field of SU ~ 1/Ca m would change |A) by adding to it a significant part of As a result, the state 
| A) is only as stable as the largest value of Ca m . The least stable eigenvector is that for which this measure of stability 
is minimized. The stability of the localization of Dirac eigenvectors in a given gauge field configuration depends on 
the least stable localized eigenvector. Hence the stability can be defined to be the quantity 



S = miiiAsioc^ 



where S\ 



(10) 



such that the minimum is over states |A) which are localized and the maximum is over states which are extended. 
The inverse, 1/S\, for a localized state |A), is a measure of the minimum field strength which causes significant mixing 
with an extended state. This measure is eminently suited to a lattice where the spectrum is discrete. If indeed there 
is stable localization, then examination of the particular element of the mixing matrix which gives S can help us to 
identify whether localization is achieved through Mott's mechanism and the formation of mobility edges, or through 
spatial segregation of the support of localized and extended states. 

A numerical implementation of eq. (|10[) requires specification of which eigenvectors are localized. We use a definition 
in terms of the IPR, taking all eigenvectors with P 2 X > P 2 * are localized and those with P 2 X < P 2 * are deemed to be 
extended. When changing this definition in the range 2 < P 2 ^ 10 wc found no significant change in the quantities 
reported below. The data shown in the figures are obtained with P 2 * = 8. 

In Figure [8l we show stability of the most localized states at 2T C as a function of P 2 . The quantity plotted is a 
dimensionless measure of the minimum change in the gauge field required to mix a given localized state with any 
extended state — 1/ATS\, As shown in the first panel, there is a tendency for 1/ATS\ at a given P% to decrease as the 
lattice size increases. Scaling the data by a power of the lattice size one finds an optimum scaling as the point where 
the Fisher's linear discriminant is least able to separate the data for different lattice sizes. In the second panel of the 
figure we exhibit the resultant scaling with the lattice size, 1/4T6>a oc (LT)~ 2 at fixed P£. If this scaling persists at 
larger lattice sizes, then it would imply that in the thermodynamic limit an arbitrarily small change in the gauge field 
can destabilize the localized eigenvalues. 

The scatter in the data does not allow us to measure the scaling exponent more precisely. One could argue that 
since we are examining localized states, the factor of pi in eq. © does not scale with volume. In that case one is 
forced to the conclusion that the observed volume dependence come from the energy differences in the denominator 
of eq. © scale as I? . Such a scaling is open to clear tests, and we perform this next. 
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Since a lattice allows only discrete eigenvalues of the Dirac operator, the origin of localization on the lattice is not a 
mystery. Nevertheless, one could try to probe the origin in more detail. In order to do this we construct two matrices 

F\» = V \/pi{r\ X)p 1 (r;n), and G A/J = tt— — , , (11) 

r l A ~Ml 

one of which, F, looks only at the spatial overlap, and the other, G, only at the overlap in energy. Using these we 
can define the notions of stability, T\ = max^exiFxn and Gx — max^ge^GA^- 

We found that S\ is strongly correlated to both T\ and Q\. As expected from the earlier argument, T\ shows no 
scaling with L. As a result, it requires no scaling when plotted against the scaled quantity L 2 /S\. On the other hand, 
whereas Q\ scales with the same exponent as S\. As a result, when plotted against the scaled quantity L 2 /S\, one 
requires the scaling L 2 /Q\ in order for the measurements to be universal. These correlations are shown in Figure [5] 
From the figures it is clear that the stability of the localization phenomenon seen at finite lattice spacing is controlled 
by the energy level differences. The situation seems to produce a curious version of Mott's argument. In this case 
we have localized states whose spatial overlap with extended states is finite. Thus there is no segregation of the 
spatial support of localized and extended states. Localization is seen at any finite volume, and is realized through the 
formation of a mobility edge. However the gap between the localized and extended eigenvalues seems to disappear as 
a power of the lattice volume. As a result there could be no localization in the thermodynamic limit. 

One can cross check this conclusion also by computing the minimum of the mobility gap; i. e., the difference between 
the maximum energy level among the localized states and the minimum energy level between the extended states. 
This mobility gap scales to zero as 1/(LT) 3 , and the scaling is not sensitive to the choice of P% used to separate 
localized and extended states in the range 2 < P 2 * < 10. 



V. CONCLUSIONS 



In this paper we have examined the eigenvalues and eigenvectors of the staggered Dirac operator evaluated on 
thermalizcd configurations obtained in simulations of QCD with two flavours of dynamical staggered quarks at tem- 
peratures between 0.75T C and 2T C with lattice spacing of a = 1 /AT. The spectrum develops a gap as one crosses T c , 
although in the high temperature phase the gap remains substantially smaller than that in free field theory. It would 
be interesting to study the gap formation in the transition region to check whether this way one can obtain additional 
insight on the crucial question of the order of the phase transition. 

The smallest eigenvalues have eigenvectors which are localized. We investigated different quantities, the inverse 
participation ratio (IPR) and the localization function, which measure the degree of localization, and found good 
agreement between them. 

We investigated the stability of localization properties of the staggered Dirac eigenvectors with respect to changes 
in the gauge field background. We showed that localization properties are not stable as one takes the thermodynamic 
limit. In fact, the scaling of the data shows that in that limit localization of staggered Dirac eigenvectors is not 
expected to be of thermodynamic importance. 

We developed measures of stability which distinguish between stability due to spatial and energy separation of the 
eigenf unctions. QCD wth staggered quarks seems to contain a curious reversal of Mott's argument. The support 
of localized wavefunctions is not spatially separated from that of extended wavefunctions, and this persists into the 
thermodynamic limit. As a result, if localization were to be obtained, it would be through the formation of a mobility 
edge. Indeed, at each volume, one does seem to observe the formation of a mobility edge. 

However, localization is spoilt by the fact that the energy denominators can become arbitrarily small, scaling as 
a power of the spatial volume in the thermodynamic limit where L — > oo. It would be interesting to extend this 
work to the overlap Dirac operator, whose exact zero modes are related to localized topological features of gauge field 
configurations [|[. 

This work was funded by the Indo-French Centre for the Promotion of Advanced Research under its project number 
3104-3. Part of the computations were carried out on the Cray XI of the Indian Lattice Gauge Theory Initiative 
(ILGTI). 
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